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Abstract: We study the possibility of obtaining large nongaussian signatures in the Cos- 
mic Microwave Background in a general class of single-field nonlocal hill-top inflation mod- 
els. We estimate the nonlinearity parameter /jvx which characterizes nongaussianity in 
such models and show that large nongaussianity is possible. For the recently proposed 
p-adic inflation model we find that /jvl ~ 120 when the string coupling is order unity. We 
show that large nongaussianity is also possible in a toy model with an action similar to 
those which arise in string field theory. 
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1. Introduction 

Though the simplest models of inflation yield a negligible degree of nongaussianity in the 
Cosmic Microwave Background (CMB) there has been considerable interest recently in 
constructing models which can give rise to large nongaussian signatures Q-Q (see Q 
for a review). Nongaussianity is typically characterized by the dimensionless nonlinearity 
parameter /tvl which is of the order of \n s — 1| <C l 1 (where n s is the spectral index) 
in conventional models Q-Q. The current observational limit is |/jvi| ~ 300 [|ll]] (for 
the WMAP 3 year data ]l2|) and future missions are expected to be able to probe Jnl 
as small as order unity |l3|]. Though it has been difficult to find models which can yield 
large nongaussianity (indeed, observation of |/jvl| ^ 5 would be strongly indicative of some 
novelty in the dynamics driving inflation) there exist several examples in the literature: 

1 Nongaussianity can further be characterized using the trispectrum, which is also small in the simplest 
models @. 
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1. 



Single field inflationary models in which the inflaton has a small sound speed [14|, as 
in (H. 



2. Single field models in which the inflaton potential has a sharp feature [16| 



3. Hybrid inflation [17] can generate significant nongaussianity during the preheating 
phase for certain values of the model parameters [18, O]. Nongaussianity from pre- 
heating has also been considered in [pC|], 

4. The curvaton model [ffif |. 

5. Ghost inflation pi. 



The simplest multi-field models do not yield large nongaussianity [23], however, this need 
not be true also of more complicated models. In this paper we show that large nongaus- 
sianity may also be generated in certain inflationary models based on nonlocal field theory. 
In particular, we will show that large nongaussianity is possible in the recently proposed 



p-adic inflation model [24]. 

Recently, cosmological applications of field theories containing an infinite number of 
derivatives have attracted considerable interest in the literature [24|-|31[. In [^] higher 
derivative modifications of gravity were considered and it was shown that such theories 
can have novel features such as improved ultraviolet (UV) behaviour and the existence of 
nonsingular bouncing solutions. In j^Tj it was shown that similar nonlocal gravity theories 
can lead to self-inflation. 

In this paper we will be interested in the cosmology of scalar field theories containing 



an infinite number of derivatives (we assume a standard gravitational sector) as in [25] 
and 



]-[31[. Nonlocal scalar field theories of this type can be derived from string theory, 
the most popular examples being the p-adic string theory 
theory f 



-[34] and cubic string field 
[37]. In both cases the scalar field is a tachyon describing the instability of 



some unstable D-brane configuration. Our interest in these types of nonlocal field theories 
is motivated by their novel cosmological behaviour. Nonlocal field theories can exhibit a 
variety of interesting cosmological properties including the possibility of slow roll inflation 

HI 25], bouncing cosmologies [29] and an equation of state w < — 1 



with steep potentials 
for dark energy within a sensible microscopic theory |§, ||, ||. 

In this paper we will focus on nonlocal field theories in which inflation is realized near 
the top of an unstable maximum of the potential (we refer to such constructions as nonlocal 
hill-top [^] inflation models). In [^4|] a particular model based on p-adic string theory, p- 
adic inflation, was proposed. A novel feature was the possibility of slow roll dynamics for 
the inflaton, despite the presence of a steep potential, which suggests an intriguing possible 
resolution for the fine tuning problems which typically plague inflationary model building. 



In [25] a more general analysis of nonlocal hill-top inflation was performed and it was shown 
that this behaviour is possible also in other models. In this paper we study the possibility 
of obtaining nongaussian perturbations in general models of the form considered in p5| . 
specializing at the end to the case of p-adic inflation. We find that in p-adic inflation 
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Inl ~ 120 is possible for g s ~ 1. We show that large nongaussianity is possible also in 
other models. 

The plan of the paper is as follows. In section ^ we consider the inflationary dynamics 
of nonlocal hill-top inflation, generalizing the results of [25| to include the possibility of a 
(f) 3 term in the scalar field potential. In section || we show how these results apply to the 
case of p-adic inflation. In section || we provide an estimate for the nonlinearity parameter 
in nonlocal hill-top inflation. In section |5] we use this result to estimate the nonlinearity 
parameter in p-adic inflation, showing that |/jvx| 3> 1 is possible. Finally, we consider 
a toy model with an exponential kinetic function (which is typical of string field theory 
Lagrangians) showing that large nongaussianity is possible. 



2. Nonlocal Hill-top Inflation 



We consider general nonlocal theories of the form [25] (see also 



C 



7 



>F 



□ 



where we assume a hill-top potential of the form 



<f>-U{<f>) 



(2.1) 



u{4>) 



(2.2) 



where the • • • denotes terms of order 0(</> 4 ) and higher. Such terms will generically be 
present, but since we are interested in the dynamics close to the false vacuum = 0, they 
will play a subleading role in our calculation. Indeed, the potential fl2.2j ) should only be 
thought of as an effective description close to = since it must be supplemented by 
additional terms to ensure that U(cp) is bounded from below — at least on one side. We 
use metric signature r]^ u = diag(— 1, +1, +1, +1) so that □ = — df + V 2 in flat space. In 
the above we take 0, Uq, fi and g to be dimensionless while 7, m s have mass dimension 
one. With no loss of generality we can set F(0) = 0, since F(0) can always be absorbed 



into the definition of fi. The differential operator F 
expansion 

□ 



F 



m- 



where 



n! 



n=0 

1 & n )F(z 



□ 
In* 

□ 

m 2 



should be understood as a series 



(2.3) 



and Co = 0. Lagrangians of the form ( |2.1| ) arise frequently in string field theory, in which 
case m s coincides with the string mass scale. In the following we keep m s unspecified. 
Examples of theories of the form (2.1) include the p-adic string theory |32j (see also 

mm) 

1 , — ^ , 1 



9 
9p 



p + 1 



if)- 



p+1 



(2.4) 
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with 



9p 



9s P 



and, in order to put fl2.4|) in the form (|2.l|), we take ■0 = + 1. In (2.4) m s and g s are 
the string mass and coupling respectively. The action (2^4) was derived assuming p to be 
a prime number, though it appears to make sense for any integer value of p. A second 
popular example is the tachyon action (see also [28, 3£, |37j ) 



ml 
C = — - 

2 

9t 



1 



□ 



mi 



1 



-0 1 + e 4m i - -0 



(2.5) 



where A 2 = 0.9556. This tachyon action is derived in cubic superstring field theory (SFT) 
in a flat background, incorporating only massless fields. 

2.1 Background Dynamics 

We first consider the homogeneous dynamics of the theory ( |2.1| ). The Klein-Gordon equa- 
tion for the scalar field is 



□ 



/?/- 



The Friedmann equation is 



0(t) = V0W + 50 2 (i) 



(2.6) 
(2.7) 



where the stress tensor is given by |25| (see also |4C|) 
1 



7' 



-T 



E 

z=i 
1 



j=0 



(d« D ^) 3, 



□ 



1 

+ 2 



i=0 



2F 



(2.8) 



In order to solve the system ( |2.6| , |2.7D we make the ansatz of a series expansion 



„r\t 



H(t) 



r=l 



r=l 



jrXt 



(2- 



and solve for A and {0 r , i^ r } order-by-order in powers of u = e xt . We have parameterized 
the solution so that as t — > — oo the field sits at the unstable maximum = and the 
universe undergoes de Sitter expansion with Hubble constant Hq. We use the freedom to 
choose the origin of time to set 0i = 1. We also set H\ = which is consistent because 
T^y does not contain any terms linear in 0. Thus, to solve the equations of motion up to 
order 0(u 2 ) the truncated expansion 



0(t) 
H(t) 



U + 2 « 

H - H 2 u 2 



(2.10) 
(2.11) 



suffices. 



- 4 - 



2.1.1 Scalar Field Evolution 

We now proceed to solve ( |2,6| ) using ( |2.10| ). It was shown in |24[ | that 

(-□)> = + (A 2 + 3H \) n u+ (4A 2 + 6H X) n <p 2 u 2 + 0(u 3 ) (2.12) 
so that the left-hand-side of ( |2.6| ) is 

F (^y F (-^y +F (-^)^ + o^ ( , 13) 

while the right-hand-side of ( |2.6D gives 

- /x 2 + g<f) 2 = -/j 2 u + (g- ii 2 4> 2 )u 2 + 0(u 3 ) (2.14) 
Matching at linear order in u implies 



It is convenient to define 

= -m 2 ^- 1 (-/i 2 ) (2.16) 

since consistency of this approach requires that F~ l exists and is single- valued, so that 
(2.15) takes the form 

A 2 + 3H X - uj 2 = (2.17) 
It is straightforward to solve ( 2.1 7| ) for A 



2 2 \/ 9^0 

The slow-roll solution corresponds to taking the positive root and assuming that uj 2 <C Hq 
so that 

A = ^rr (2.18) 



We also define a dimensionless parameter 

,2 



77 = — ^ (2.19) 
3#o 



so that slow roll corresponds to |??| <C 1 (notice that r\ < in this model). For the slowly 
rolling solution we have X/Hq = \r]\ <C 1 and A 2 <C 3f^oA so that the evolution is friction- 
dominated in the usual sense, namely <p <C H<p when \rj\ <C 1. In passing, notice that to 
leading order in the small-u expansion the field obeys the eigenvalue equation 

□ = -uj 2 (j) (2.20) 

However, the correspondence between the solutions of ( p. 20 ) and the solutions of (| 
breaks down beyond leading order in the small-u expansion. 
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We now solve (|2.q ) up to second order in the small- w expansion. Notice that the 
arguments of F in ( |2.13| ) are approximately —ui 2 /m 2 and —2uj 2 /m 2 for the coefficients of 
the 0(u) and 0{u 2 ) terms respectively. Matching ( [2.13 ) and ( 2.14 ) at order u 2 gives 



H 2 + F 



(2.21) 



for \rj\ <C 1. Notice that in the case g = we have (f>2 = and the correspondence between 
the solutions of ( 2.20|) and the solutions of fl2.6|) holds quite generally. As we have pointed 
out previously, however, for g ^ this equivalence breaks down at order u 2 . 

2.1.2 Friedmann Equation 

It now remains to determine the coefficients Hq, H 2 in ( [2.11 ). In order to determine H(t) 
up to order u 2 we need only consider <j)(t) up to first order in the small-u expansion. This 
is so because T^ u does not contain any term linear in </> (see eq. fl2.8|) ). Working strictly 
with 4>{t) = u the scalar field obeys ( 2.20| ) so that it is straightforward to resum the infinite 
series in (|2.8|). The result is 



-4 rp 



" 2 F' (-£ 

m 2 



2m 2 



i_ H 
3 



u 2 + 0(u 3 ) 



Then we have 



M)0 



V3M p ^M p 



4m 2 U ~ 



mi 



F'l-^)u 2 + 0{u*) 



(2.22) 



(2.23) 



at leading order in the 77 slow roll parameter. The fact that there is no 0{u) term on the 
right-hand-side of ( [2.22 ) demonstrates that it was consistent to set H\ = in ( p. 11 ). From 
the Friedmann equation 

Hit) = -^«L 
V3M P 



it is straightforward to identify H , H 2 by equating (2.11) to (2.23). We find that 



y/SM v 



and 



th _ 

H AmlUv 



F' 



mi 



Or, using ( 2.16 ) to eliminate the derived parameter u , we have 



t = -4^ F - I (V)F'[F- 1 (V)] 



H, 

It is useful to define a second dimensionless "slow roll" parameter by 

Ho 



e(t) 



2Xt 

H 6 



(2.24) 



(2.25) 



(2.26) 



(2.27) 
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It is clear that once e(t) = 1 then the perturbative expansion ( |2,11| ) breaks down and our 
solution is no longer valid. At this point one may assume that inflation has ended and 
hence the condition e(t en d) = 1 implicitly defines the time t en d at which inflation ends. It 
is worth pointing out that the definition ( 2.27] ) differs from the usual e slow roll parameter 



since 



H 
IP 



2\r]\e 



with our definitions ( |2 .19 , [2.27 ). A quantity in which we will be interested in later on is 
the value of u = e A * at the end of inflation. Using ( |2.26| ) we see that 

1 1 



F-\-^)F' [F-\-^)] 



(2.28) 



2.2 Fluctuations 

We now study the scalar field fluctuations about a homogeneous background, 

0(t,x) = 0o(t) + 50(t,x) 



(2.29) 



As we have shown the background solution during inflation is described by 4>o(t) = e xt , 
H(t) = Hq — fye ■ In order to solve for the fluctuation 5(f) we treat the background 
as pure de Sitter space (which is equivalent to working to zeroth order in the small-n 
expansion so that (j)Q = and H = Hq). In this limit the fluctuation 5(j) obeys the equation 

F (^] 5<t> = -ii 2 5(t) (2.30) 

We can obtain solutions of ( 2.3C| ) by choosing 8<j) to satisfy the eigenvalue equation 



-oj 2 5(j) 



(2.31) 



with uj 2 given by (|2.16| ) as above. The solutions of (2.31) are well known. However, to 
make contact with the usual treatment of cosmological perturbations we need to define a 
field in terms of which the action looks canonical. This presents a difficulty because, in 
general, there is no local field redefinition which will bring the kinetic term (j)F ( ) <f> into 



canonical form. However, as in p4|, |2q], we can circumvent this difficulty by noticing that 



(2.30) can be derived from the perturbed Lagrangian 



XL 

2 



54>F 



□ 



(2.32) 



Then, for on-shell fields (that is, when ( |2.31| ) is satisfied), 



C 



(2) 

on— shell 



60F I ~ 



rat 



1 



(-u 2 )' 
1 

(-0J 2 ) 
UJ 2 



r_ 

2 

± 
2 

2 

-6<pD5<p + y(M 2 



m 



2 



111- 



□ 



-Lu 2 )5(p + fi 2 {5<j)) 2 
50 + /x 2 (^) 2 
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where we have defined 

Sip = A5<p (2.33) 

with 

(— uj l ) \ nig J rUgi* {— \x ) 



In the second equality in fl2.34|) we have used ( 2.16 ) to eliminate the derived parameter uj 



Thus Sip is the variable which has canonical kinetic term in the action. 2 

Notice that the definition of the canonical field, eq. ( 2.33j ), has been derived by study- 



ing the linearized theory. However, we will continue to adopt this definition even up to 
second order in cosmological perturbation theory when we compute the nongaussianity in 
section ||. This is justified since in the ADM formalism |3J| it suffices to use the free- 
field solution to compute the interaction Hamiltonian. This is so because the terms in 
the Lagrangian which would provide nonlinear corrections to the free-field dynamics are 
always multiplied by terms proportional to first order equations of motion, and hence these 
corrections vanish [111] . Similar comments apply to our construction of the curvature per- 



turbation in subsection [2.3[ In any case, contributions to /jvx coming from the free- field 
dynamics must always be present and, barring anomalous cancellation, the contribution of 
such terms provides a lower bound on the actual nongaussianities produced. 



We now proceed to solve ( 2.31 ), bearing in mind that Sep is the appropriate canonically 



normalized field. We write the Fourier transform of the inflaton fluctuation as 

8y(t,x)= [ _^!*- e *-£ k (t) (2.35) 



(2tt) 3 /2 

where the operator- valued Fourier coefficients can be decomposed into annihila- 

tion/creation operators a/^aj, and c-number-valued mode functions tp^ as 

= a k p k {t) + at- k ¥>-k(t) ( 2 -36) 
The mode functions <pk(t) are given by 

(2 - 37) 

where the order of the Hankel functions is 

[9 ^~ 

v =\ 7 + ^53 ( 2 -38) 



and of course a = e Hot . In writing ( |2.37| ) we have used the usual Bunch-Davies vacuum 
normalization so that on small scales, k 3> ciHq, one has 

,-1 



a 



2k 



2 Our choice of normalization coincides with the definition of a "canonical" inflaton advocated in [ pi| . 
However, this choice differs from the normalization employed in |25|, where the inflaton was normalized in 
such a way that the stress tensor takes canonical form, though the kinetic term in the action does not. 
In the case of p-adic inflation, our main interest, the discrepancy is only a factor of \/hxp which is of order 
unity for the values of p which we consider. 
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which reproduces the standard Minkowski space fluctuations. This is the usual procedure 
in cosmological perturbation theory. However, we note that the quantization of the theory 
fl2.1| ) is not transparent and it might turn out that the usual prescription is incorrect in 
the present context. We defer this and other subtleties to future investigation. On large 



scales, k <C alio, the solutions (2.37) behave as 



H 



\tPk\ 



/2P 



k \ 3 / 2 ^ 



which gives a large-scale power spectrum for the fluctuations 

2 / , \ n s -l 



Ho 

2tt 



k 



aH, 



with spectral index 



7U 



l = 3-2u = 2r) 



(2.39) 



where rj is given by ( 2.19| ). Since rj < this model always gives a red-tilted spectrum, in 
agreement with the latest WMAP data Q . 

Using (2.19) and ( 2.24| ) we can rewrite ( |2.39| ) as 



Us 



Uo 



\F~ 



(2.40) 



which can be use to eliminate 7 in favour of other parameters. 

It is worth noting that the equivalence between the solutions of ( 2.3 1| ) and ( 2.3C| ), 



which implies the equivalence between local and nonlocal theories described in [25], appears 
only at linear order in perturbation theory. Beyond linear order (as long as g 7^ 0) this 
equivalence breaks down, as we illustrate for the case of p-adic inflation in appendix A. The 
breakdown of this equivalence, once interactions are included, is crucial for understanding 



why nonlocal theories of the type (2.1) can give rise to significant nongaussianity. Indeed, 
if it were true that the equivalence persisted at all orders in perturbation theory then the 



calculation of /nl in the theory (2.1) would be exactly equivalent to a calculation of f^L 
in some local field theory where f^L <C 1 is generic. 

2.3 The Curvature Perturbation 

A full calculation of the power spectrum should include metric perturbations and also de- 
viation of the background expansion from pure de Sitter. Such a computation is beyond 
the scope of the present paper. Here we neglect metric fluctuations and take H = Hq to 
compute the field perturbation, as we have done in the last subsection. (These approxima- 
tions reproduce the full calculation up to acceptable accuracy in the standard theory.) We 
assume that the curvature perturbation is given by 

H 



c 



-S(p 



(2.41) 



where (po = A<po = Au. To evaluate the prefactor H/ipo we must work beyond zeroth order 
in the small-it expansion. We take <po = u to compute H/(fio even though the perturbation 
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5<p is computed in the limit that 0o = 0. This should reproduce the full answer up to 0(u) 
corrections. It will come in handy later on to define £ = c^5ip so that the prefactor is 

c =-— = -^- 

tfo AX u 

To compute the spectrum this must be evaluated at horizon crossing, N e e-foldings before 
the end of inflation. Noting that 

u(t) = a(t)M 

and the scale factor at horizon crossing is related to that at the end of inflation by a+ = 
e~ Ne a en d, it follows that 

u* = e' N ^u cnd = e-tl"--\ cnd (2.42) 

while u en d is given by ( p. 28 ). In the above = u(t+) = e xt * is the value of u(t) at the 
time of horizon crossing, t*. Further, since \/Hq = \rj\ = \n s — l|/2 and using ( |2.34| , |2.28 
2.42j ) we can write the prefactor = —H/ipo at the time of horizon crossing as 



2 _( H\ 2 4 rr^F- 1 (-n 2 )e N ^ n '- 1 \ 

C/- 



c \ K-l| 2 7 4 (-/^ « 2 



fi- 



end 



N e \n s -1\ 2 -1 

2^fr\ [F-H-V 2 )] 2 F' [F-H-» 2 )] (2.43) 



\n s - l\ 2 7 4 C/ \i 

it follows that the power spectrum of the curvature perturbation is 

n,- 1 



where 

^K-IQW^w 325 * 10 ' 10 ^ 

Using ( 2.43| ) and ( |2.24| ) we arrive at the result 

where we have defined the function 

G{f) = \ [F-\-f)fF> [F-\-^)] (2.47) 

3. Specialization to p-adic Inflation 

We now consider how our previous results apply to the case of p-adic inflation fl2,4| ) . For 
the dimensionful parameters we have 

4 4 2 

9p gj P - 1 
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and m s is identified with the string mass scale. Notice that 7 > m s for typical values of g s 



and p. The dimensionless coefficients appearing in the potential (2.2) are 

lp-l 



U 
M 2 



2p+ 1 
p-1 



P, 

- 2 {v 



f'l --(/'-I) 

The function F associated with the kinetic operator is 



F(z) = 1 - p 



-2/2 



so that 



F' [F-\- 



P 



Inp 



G(fi 2 



2 

2plnp 



P 



1 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



(3.6) 
(3.7) 

(3.8) 



Notice also that ( 2.4C| ) can be used to eliminate the ratio m s /M p in favour of p, g s as 

8(p + 1) 



m s 

M r 



St 



Using equations (pT6 -3.8) it is easy to show that we exactly reproduce the expressions 
for A, fa, Ho, H2 and n s — 1 given in |f24| |. However, our results for u en d and A 2 (equation 
( 2.28 ) and ( 2.46[) ) differ somewhat from the results of [24]. To see this we first compute 



Uend for p-adic inflation using (|2.28[) 



A end 



l p+ 1 

2p- 1 



lnp 



Now, using ( |2.46| ) to compute the COBE normalization for p-adic inflation we find 

lp+1 



A- 



8 g;e 



2 „JV e |n s -l| 



37T 2 p 



II 3 



2p- 1 



Inp 



(3. 



(3.10) 



Equations ( |3.9| ) and ( 3.10| ) should be compared to equations (5.18) and (5.20) in 24]. In 
both cases the equations differ by a factor of (p+ 1) Inp/ [2{p — 1)]. This discrepancy arises 
because the end of inflation was defined here as the end of the validity of the perturbative 
expansion is powers of u, whereas in [24] the end of inflation was defined using the friction- 
dominated approximation, which has a longer range of validity. (This discrepancy was 
noted also in ]|25|].) Thus, in the particular case of p-adic inflation, the slow-roll dynamics 
actually persist even after the perturbative description ( |2.9| ) has broken down. Thus, 
choosing ( 2.28| ) to define the end of inflation is more stringent. 

Notice that in the case of p-adic inflation this discrepancy does not lead to any signifi- 
cant quantitative change in the results since the factor (p + 1) Inp/ [2(p — 1)] is order unity 
for the values of p which were considered in [24], and which will also interest us here. 
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4. Calculation of the Nongaussianity 



We now proceed to estimate the level of nongaussianity in nonlocal models, based on the 
bispectrum. As explained above, a full calculation including metric contributions appears 
to be prohibitively difficult; we do not know how to quantize the metric perturbations due 
to the complexity of the stress-energy tensor. However using local scalar field theories as a 
guide, we hypothesize that the scalar contribution to the curvature perturbation is of the 
same order of magnitude as that of the metric contribution, so that a valid estimate can 
be obtained from the scalar perturbation alone. Its contribution to the bispectrum comes 
from the cubic term in the action for small fluctuations. 

4.1 Perturbing the Lagrangian 



It is straightforward to perturb the Lagrangian (2.1) up to cubic order in the field Sip = Adcj). 
The quadratic Lagrangian is 

,2 



UJ 

V 



8<pF[ ^)6<p + fJ ?(6ip) 2 



The cubic Lagrangian is 
where we have defined 



= c H {5<pf 



CH 



gm° s 

3 7 2 



2^3/2 



(4.1) 

(4.2) 
(4.3) 



4.2 Calculation of the Bispectrum 



Variation of the quadratic Lagrangian ( [4.1| ) yields the free perturbation equation of motion 
(2.30) whose solutions are given by ( pBl p6| , |l37l ), as we have already discussed. To 
leading order in slow-roll parameters the order of the Hankel functions in ( 2.37| ) is v = 3/2 
so that 



¥>k(r) 



~ AS H 



-ikr 



(1 + ifcr) 



(4.4) 



where we have employed conformal time r, which is related to cosmic time t by adr = dt. 
It follows that a = e Hot = — 1/(Hqt). In ( [4.4D k is the comoving wave number of the 
perturbation and 5 is a constant real phase whose value is irrelevant for our calculation. 

The quantity of interest for our calculation is the 3-point correlation function of the 
curvature perturbation (CkiCk 2 Ck 3 ) where the expectation value is computed in the inter- 



action vacuum e dt Hint ^ ^\0). The leading contribution is 



(C kl Ck 2 Ck 3 (r)) = -i / dr'(0\ [Ck 1 Ck 2 Ck3(r),^int(T , )] |0) 

= -icg [ T dT'(0\ [^ 2 ^ 3 (r),H iQt (T')] |0) 
J n 

where we have used the fact that 

Ck(r) = c ? ^ k (r) = c c a k v3 k (r) + aL k V-k( T ) 



(4.5) 
(4.6) 
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The interaction Hamiltonian takes the form 

H int = ~ ! d 3 xa 3 C^ (4.7) 



c H I d 3 xa 3 5(p 3 (4.8) 



where ch is defined in ( [4.31) . 

In order to compute (fO|) we first consider the commutator 



g— ik-x 



[e k (r),^(r',x')] = —r^ (^(t)^(t') - ^(r')^(r)) 



ff2 -ik-x' 



3 v y (2vr) 3 / 2 

where on the last line we have taken the large scale limit —kr, —kr' <C 1. 
Noting that the commutator (fOI) is a c-number we have 



[6c(t),«V(t',x')] = 3 [£ k (T),<V(r',x')] W,x') 

The expectation value in (^1]) can be written as 

<0| [& 1 &MT),8ip 3 (r',x!)] |0> = 3^ k3 (r),^(r',x')] (OI&^t) W, x') 

+ 3[£ k2 (r),M<x')] <0|£ kl <^V,x')£k 3 (r)|0) 
+ 3[£ kl (T),ty>(T' J x')] (0|V(^x%^(t)|0>(4.10) 

Consider, for example, the first term. Carrying out the Wick contractions we have 

3 [e k3 (r),^(r',x')] <0|£ kl £k 2 (T)^V,x')|0> 
= 6 [e k3 (r),^(r',x')] e'^+^D^r - t')D^{t - r') 

where the momentum space propagator is 

1 H 2 

e- ik{T - T '\l + ikr)(l - ikr') 



(2vr) 3 /2 2k 3 
1 H 2 
(2vr) 3 / 2 2P 



(4.11) 



and on the last line we have restricted to large scales. The remaining terms in ( [4.10D are 
similar. 

We are now in a position to calculate ([4.5[) . On large scales — k{T, —kiT 1 < 1 we have 
(C^OA (r)) ^ -c H cl^ £ dr> (T " T } '/ k * + | + k * 5® (kl + k 2 + k 3 ) (4.12) 
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The time integral is divergent and needs to be regulated. Taking r to correspond to the 
end of inflation and t% to correspond to the beginning of inflation we have 

dr' (T ~ T ? 3 In ( T S\ = In e N « = N e (4.13) 



(t') 4 \t 

where N e is the total number of e-foldings of inflation. 

It is worth commenting on the infrared (IR) divergence in the dr integral in ( 4.12| ) which 
does not appear in the standard calculation of Jnl H> @- Naively, this divergence would 
seem to suggest that the curvature perturbation in £ is not freezing out on super-horizon 
scales, as it does in the standard theory. However, we suspect that the IR divergence in 
Q4.13 ) is in fact an artifact of the approximations which we have made. Indeed, a calculation 
of Jml in a local field theory which makes precisely the same approximations as we have 
made (neglecting the metric perturbations and the departure of the background from pure 
de Sitter space) yields exactly the same IR divergence |^]J . It is clear from the calculation 
of |44]] that this IR divergence will generically arise in the three-point function of any light 
scalar field in de Sitter space. In appendix B we show that, despite the presence of the IR 
divergence, the calculation of [41| can be used to estimate /nl in a local field theory and 
their result does reproduce the result of the complete calculation || up to factors of order 
unity. We assume that the same holds true in the nonlocal theories which we consider. 

Though we believe that in a more complete calculation one would find that £ is con- 
served outside the horizon in the nonlocal theory, we were not able to demonstrate this 
conclusively. A proper resolution of this issue would require us to perform a complete 
second order cosmological perturbation calculation, including metric perturbations and 
departures of the background from pure de Sitter space. Such a calculation is compli- 
cated by the nonlocal structure of the theory and is beyond the scope of the present work. 
Though it is a logical possibility that £ is actually not conserved in the nonlocal theory, 
it would be surprising. Quite generally, the time derivative of the large scale curvature 
perturbation is proportional to the nonadiabatic pressure, £ cx P na d [42]. In a model with 
only one degree of freedom there are no entropy perturbations and P na d = so that £ = 
on large scales. Thus we expect that as long as (2A) describes only a single degree of 
freedom then the large scale curvature perturbation should be conserved (in a complete 
calculation). It is well known that higher derivative theories with finitely many derivatives 
can be rewritten in terms of multiple scalar fields (at least one of which will generically 
be ghost-like). However, for infinite derivative theories of the type which we consider this 
need not be the case j^6j. For infinite derivative theories one should consider the pole 
structure of the propagator Gk ~ \F{—k 2 /m 2 s ) + /i 2 ] . Our restriction to models where 
-F _1 (— /x 2 ) is single- valued (below equation (|2.16|) ) ensures that the propagator has only 
one pole (at least near the false vacuum 4> = 0). 3 It follows that our analysis is restricted to 
infinite derivative theories which describe only a single scalar degree of freedom and hence 
one expects _P na d = (and therefore £ = 0) for such theories. 

For example, in the case of p-adic inflation the propagator near the unstable maximum has only one 
pole corresponding to a scalar field with mass-squared —2m 2 s : the bosonic open string tachyon. In contrast 
there are no poles about the true minimum, which is the p-adic version of the statement that there are no 
open strings in the tachyon vacuum. 
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4.3 The Nonlinearity Parameter 

We now estimate the nonlinearity parameter by comparing ( [4.12| ) to the WMAP ansatz 
(C kl Ck 2 C k3 (r)) = (2vr) 7 (-^/jvl) ^l^ 1 ^ (kx + k 2 + k 3 ) (4.14) 

where, as we have shown previously, A 2 = H 4 /(2-Kipo) 2 = c 2 H^ / '(47r 2 ) . Comparing ( 4.14| ) 
to ( [4.12| ) we have 

5iV e 1 

fNL = 3p W ^ (415) 
It is straightforward to compute the product Cf/c^ using (f4.3|) and ( 2.43| ). The result is 



3Uq /2 V * J I"- - M V [ -l>- 



where G(/j?) is defined in ( 2.47|) . It is possible to further simplify this expression by using 
(2.40) to eliminate 7 -4 with the result 

^ = -g^feV^K- 1 |g^=g!>^ygP) (4 . 17 ) 



6 \M P J |//| 3 



We can eliminate (m s /M p ) 2 from ( 4.17 ) using the COBE normalization ( 2.46 ). The final 
result is 

fNL = -^gUi' 2 e-**>-%. - 1| 2 K(^) (4.18) 

where we have defined 

K{p2) = i^nyyp^ [i ?_i (-^ 2 )] i 1 / 2 (419) 

Equation ( (4.18|) is the main result of this section. 
4.4 The Perturbative Regime 

Since fNLOtg and we are interested in models with |/jvl| ^ 1 it is important to understand 
what values of g are acceptable. In conventional (local) inflationary models one would 
require \g\ <C 1 in order not to spoil slow roll; however, this need not be the case for the 
nonlocal hill-top models which we consider. Indeed, in the case of p-adic inflation one has 
\g\ ~ p 2 S> 1 (see section |3|). In terms of the canonical field the Lagrangian is of the form 

£ = ^ F (4)"- FW <420) 

where 

V(ip) = V -^ 2 + c H( p 3 (4.21) 

and Vo = 7 4 f7o- In a local field theory the one-loop correction to the cubic term of the 
potential ( [4.21 ) would be of the order 

16tt 2 lo 3 
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so that the dimensionless quantity Cjj /u> 2 controls the perturbative expansion. We therefore 
expect that the calculation is under control as long as c 2 H < uj 2 . Using equations ( [Q| ) and 



( 2.16 ) the condition (? H < u implies the following upper bound on g: 



\9\ < 



3/i d 



1_ 
m.. 



(4.22) 



In a given model the largest reliable value of Jnl is achieved by choosing g to saturate the 
bound flOij ). 



5. Examples 

5.1 p-adic Inflation 

With the results of section ^ we can compute /jyx for the p-adic inflation model, using 
pTTSj p^D. The result is 



/iVL 



5iV e 
1927T 3 



r e 2 



|n s -l] 



n.. 



l| 2 v^ 



2 p-1 
lnpp + 1 



1/2 



so that Jnl ~ a/P for p » 1. This result may slightly underestimate the nongaussianity 
produced in p-adic inflation because of the extra factor 2{p—\)/[(p + 1) lnp], which appears 
because the perturbative approach breaks down before the apparent end of the slow roll 
phase. This is precisely the discrepancy discussed below ( |3.10 ), Taking into account the 
fact that in p-adic inflation slow roll ends when u ~ p^ 1 ^ 2 , rather than when u is given by 
(|3.9| ), we obtain 



NL 



5N t 
192^r 3 

For n s ^ 0.95 and N e ^ 60 this implies 



f NL 2.8 x 1(T 5 



(5.1) 



(5.2) 



so that large f^L requires large values of p. 

The CO 

ilO^e-^l^-^ln 



How large can (|5.2| ) be made? The COBE normalization relates g s and p as [24| 

5tt\/3 



9* 



2\/2 



l| 3/2 VP 



For N e = 60 and n s = 0.05 the condition g s < 1 bounds p from above as 

p < 1.7 x 10 13 



so that, using (5.2), the nonlinearity parameter is bounded as 



fNL £ 120 

with the upper limit corresponding to g s = 1. The maximum possible nongaussianity is 
within the observational limit fNL < 300 though it should be observable in future missions. 
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We see that large nongaussianity is possible in p-adic inflation. It is interesting that the 
largest values of /jvx correspond to g s close to unity, which is considered natural from the 
string theory perspective. 

In evaluating out estimate for Jnl, eq. (|^), we have used the total number of e- 
foldings between horizon crossing and the end of inflation N e = 60, rather than the number 
of e-foldings which can be observed. This is justified because the factor N e has its origins 
in an IR divergence ( pj~l3| ). This diver gence should be regulated by the largest scale to 



which we have experimental access, namely, our current horizon size. The need for such IR 
regulators in cosmological perturbation theory was discussed in Q where it was argued that 
these divergences simply reflect our ignorance about scales beyond the horizon. Notice that 
typically the total number of e-foldings between horizon crossing and the end of inflation 
is iV e ~ \n s — so that one of the factors of n s — 1 in ( |4.18| ) is cancelled by the factor 
N e . An identical cancellation occurs in local field theory models, see appendix B. 

Finally, we note that even taking N e = 6 we would still obtain f^L ~ 12#s which, 
for g s ~ 1, is still several orders of magnitude larger than the prediction of the simplest 
inflationary models. As long as Jnl > then the primordial nongaussianity should 

be detectable in future missions p3| . For the primordial /jvl = 0(1) the situation is 
more complicated since post-inflationary super-horizon evolution will also generate an order 
unity contribution to the nonlinearity parameter JIJ (see also [|j ) and it may be difficult to 
disentangle these two distinct sources of nongaussianity. 

In principle one could make Jml even larger than /jvz, ~ 120 in this model by relaxing 
the requirement that g s < 1, however, at large g s we can no longer rely upon perturbation 
theory to compute the bispectrum or other quantities from the tree-level p-adic action. 
String loop effects will invalidate this action when g s > 1. This is also consistent with 
our discussion in subsection 4.4 since the bound CH), for the case of p-adic inflation, 
corresponds to g s < 3. 

5.2 A Toy Model with Exponential Kinetic Function 

As an example of further possibilities within the framework of nonlocal field theories, we 
next investigate a toy model which can give rise to significant nongaussianity. We choose 
a kinetic function 

F(z) = - (l - e-f") (5.3) 



a 

where it is assumed that a,/3 > 0. Equation ( [slf ) is typical of the kinds of kinetic functions 
which arise in string field theory. It is straightforward to compute 

F~\-» 2 ) = -^ln(l + a^ 2 ) (5.4) 

F'[F-\-f)] =^(l + a/x 2 ) (5.5) 
a 

so that 
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1 + aji 2 y In (1 + afj, 2 ) 



(5.7) 



In passing, it is interesting to compute the effective mass of the field fluctuation for this 
model using (|[~J) and ( |2.16| ) 



The fluctuations Sip behave as though the field had mass-squared —uj 2 . However, for an 
exponential kinetic function (|5~J) this effective mass is only logarithmically sensitive to the 
actual mass \x. This explains the novel behaviour first noted in p3] that the cosmology of 
the p-adic tachyon is virtually insensitive to the naive mass of the field. 
We set 



m 



7 



9. 



(5.8) 



by analogy with the D-brane tension. Taking Uq = 1 the equation (2.40) fixes the string 
scale as 



M r , 



while the COBE normalization requires 



2g 2 ln(l + 



\n, - 1 



P 



1 



4 = 6tt 2 K - ife-^^-^r—, 

(3 2 11 1 + a^ 2 In 3 (1 + a/i 2 ) 



• 25 • 10 



-10 



If we take a = jj, = 1 then (5.10) requires 



^ = 10^ 



(5.9) 



(5.10) 



(5.11) 



and ( |5.9| ) gives (m s /M p ) 2 = 3 x 10~ 5 g s so that m s <C M p unless g s is very large. 

To obtain an upper bound on f^L we note that g is bounded from above by ( 4.22j ) so 

that 

3fi 3 (3 



9 < 



3 ^ = 0.4x 10 7 



ln(l + a[i 2 )g s ln2 g s 
(we are taking g > 0). Using flO]]), ( pU2D , n s = 0.95 and iV e = 60 we obtain 

I/tvlI =9.5 x 10~ 5 5 



(5.12) 



(5.13) 



so that /jvx ~ 300 bounds 5 from above as g <, 0.31 x 10 7 . Since this upper bound is 
stronger than the bound ( |5.12 ) it follows that it is possible to saturate the observational 
limit on f NL . 



6. Conclusions 

In this paper we have studied the nongaussianity produced during nonlocal hill-top infla- 
tion. We were particularly interested in p-adic inflation, in which case large nongaussianity 
is possible, and indeed is natural since the upper bound on /jvx corresponds to g s ~ 1. 
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We also considered a toy model with an exponential kinetic function, typical of string field 
theory Lagrangians, and showed that for certain values of the parameters large nongaus- 
sianity is possible in this model. Thus, nonlocal hill-top inflation models are among the 
few inflationary scenarios which can give rise to /jvx S> 1. 

There are several caveats to our work. In computing the curvature perturbation £ 
we have neglected both metric perturbations and also departures of the background from 
pure de Sitter space. Although these approximations reproduce the correct answer up to 
acceptable accuracy in standard (local) theories, it is not clear if this is also true in the 
nonlocal theories which we have considered. A cause for concern is the IR divergence ( 4.1)j ) 
which we have regulated by the number of e-foldings of inflation. This divergence does not 
occur in the standard theory and seems to indicate that our definition of £ does not freeze 
out on large scales, as it should. We speculated that this divergence is an artifact of the 
approximations which we made, and this claim is supported by the fact that an identical 
divergence occurs also in the local theory when one makes exactly the same approximations 
which we have made [41]. That being said, we have not proven that £ will really freeze out 
in a more comprehensive treatment. It is possible that the relation ( ~ —HS(p/(po does not 
hold for nonlocal theories. In order to answer this question in a satisfactory manner it will 
be necessary to perform a complete and rigorous analysis of the cosmological perturbations 
in this theory, incorporating also metric perturbations. Such an analysis is complicated 
by the fact that the kinetic operator F(0) will be of the form F(O + So) where O 
is the covariant d'Alembertian associated with the homogeneous background and So is 
the perturbation which contains both metric perturbations Sg^ u and also derivatives. In 
general O and SO do not commute, making it extremely difficult to solve the nonlocal 
evolution equations. We leave a rigorous treatment of the cosmological perturbations to 
future analysis. 

In a local field theory the gravitational sector is known to give a small contribution 
to the nongaussianity and hence the neglect of the metric perturbations is justified when 
Inl 3> 1 101 • ^ * s n °t c ^ ear ^ t ms is also the case in nonlocal theories. In a complete 
calculation one would find terms in the perturbed Lagrangian which involve the operator 
F(O /m^) acting on the scalar metric perturbations and hence it is not clear that the 
contribution of such terms to Jnl will be suppressed. Barring anomalous cancellations we 
would expect such terms to be of the same order of magnitude as those which we consider 
and hence our estimate should be reliable up to factors of order unity. 

Another caveat is the complication in setting up the initial value problem for differential 
equations with infinitely many derivatives. It is known that the initial value problem for 
equations with infinitely many derivatives is fundamentally different from the initial value 
problem for an equation with N derivatives where N 3> 1 pq ] . Infinite derivative equations 
have been studied from a mathematical physics perspective in [43] . For further discussions 
about the initial value problem see p4j, H3] . 

A related difficultly is the issue of quantizing a nonlocal theory. In our analysis we have 
assumed the usual Bunch-Davies vacuum for the perturbation Sip, though it is not clear that 
this is correct in the present context. Quantization of nonlocal theories with finitely many 
derivatives generically leads to ghost excitations (and also to a classical pathology called the 
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Ostrogradski instability) though nonlocal theories with infinitely many derivatives (such 



as we consider) can evade these difficulties. For further discussion see [24, 27]. 

Despite these difficulties, we believe that our calculation does provide a rough (order 
of magnitude) estimate of the actual nongaussianity produced. It is easy to intuitively 



understand why /jvx can be made large in these kinds of models. Our result ( 4.1§| ) is 
equivalent to computing the three-point function for a light scalar field in de Sitter space. 
This correlator is proportional to coefficient of the cubic term in the potential, g, and hence 
will be large when g is large. In conventional inflationary models this is impossible since a 
large coupling g would spoil slow roll. However, in nonlocal theories this is not necessarily 
true. (Similarly a large mass term fi 2 would spoil inflation in a local theory, though there 
is no problem with exponentially large fj? in our models.) Indeed, in the examples where 
we have found large nongaussianity this is simply because g > 1, and the novelty is that 
such theories can support inflation. 
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APPENDIX A: p-adic Scalar Field Evolution Beyond Linear Order 

In this appendix we demonstrate that the equivalence between local and nonlocal theories 
breaks down beyond linear order in perturbation theory by studying the p-adic scalar field 
equation of motion 

p- n / 2 ^j = ^ (A-l) 

(in units m s = 1, which we employ throughuot this appendix) up to second order in 
perturbation theory. For simplicity we neglect metric perturbations and assume pure de 
Sitter background. Expanding the p-adic scalar in perturbation theory as 

i/> = 1 + 

= 1 + 5(1)0+1^)0 

It is straightforward perturb the field equation ( |A-1| ) up to second order with the result 

5^<f> = (A-2) 
5^4> = p(p-l) (S (1) 4>) 2 (A-3) 



-□/2 
p 1 — p 



-□/2 
p 1 — p 



At linear order one may construct solutions of ( A-2j ) by taking OS^<f> = —uj 2 5^4>, however, 



this prescription fails for the second order equation (|A-3| ) . It is interesting to notice that 
in terms of the canonical field 
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the second order equation ( A-3 ) becomes 



<5( 2 V = c(<5 (1 V) (A-4) 



where 



—n/2 
p 1 — p 



The COBE normalization gives g s ~ 10 \fp (for n s = 0.95 and N e = 60) so that c is of 
the order 

c ~ 10- V /2 (A-5) 

for p> 1. We see that the second order effects can be made large by taking p> 1 which 
is precisely the regime in which /tvl > 1 occurs. As one might expect, this also coincides 
with the regime in which the nonlocal structure of the theory is playing an important role 



in the dynamics, as emphasized in 24]. 



APPENDIX B: Comparison to Local Theory 



In [41] Falk et al. studied the three-point function of the inflaton perturbation making 
the same approximations as we have made in our calculation. Namely, |^l[ neglect metric 
perturbations and departures of the background expansion from pure de Sitter. For the 
potential 

V(<p) = V - (B-l) 

Falk et al. find 

?tt 3 V 3 k 3 ( 3 \ 

(^W^hW^W) =+— gH 2 N e ^^6^ P>> (B-2) 

i=l i \i=l / 

The factor of N e in ( |B-2[) arises from regulating a time dependence of the form lnr, exactly 
as in ( |4.13 ). We now use this result to estimate /nl for the potential ( |B-1| ), showing that 



the answer agrees with a more careful calculation || up to factors of order unity. For 
simplicity we assume that g > throughout this appendix and restrict ourselves to the 
case where ip rolls from the unstable point <p = towards some positive value tp > 0. 

We assume that Vq S> gip 3 throughout inflation so that H = Hq = \/Vq/ 'y/3M p . The 
slow roll parameters, e = 2" 1 Mp(V' /V) 2 and r\ = M 2 V" /V , evaluated at the time of 
horizon crossing t = t*, are 



p v " m 2 

3 

where (p± = <p(t = t+) is the value of the inflaton at horizon crossing. The spectral tilt is 

n s - 1 = 2?? - 6e =i 2r] (B-3) 
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The slow roll Klein-Gordon equation 



3H(j> 2i-V' 



(B-4) 



has solution 



<p(N) ^ 



1 - i-liN 
2 



(B-5) 



where iV = Ho(t — £*) so that <p(N = 0) = is the inflaton value at horizon crossing. It 
is straightforward to see that 



1 <P 2 ^ 



2M p 2 i/ 2 8M| 



so that inflation ends iV e e- foldings after horizon crossing when \rj\N e /2 ^ 1. It follows 
that \n s - l\N e = 0(1). 
We assume that 

C = c^Sip = — r 5ip 



9 



The value of cq at horzion crossing is 



H 





N-- 



=0 c^/diV 



7V=0 W* 



where in the last equality we have used flB-5|) . The COBE normalization is 



(B-6) 



(B-7) 



so that, using ([B-7D , we have 



A\ -«?(§) =25X10-" 



4 2 ~ 

C -2^2, n 2 



(B-* 



(B-9) 



Using (B-6) and ( p3-7| ) we can convert (B-2) into an estimate for the bispectrum of the 
curvature perturbation 



(B-10) 



vi=l 



Comparing (B-10) to the WMAP ansatz (4.14) we obtain the following estimate for the 
nonlinearity parameter 



\f NL \~^(n s -l) 2 N e 



(B-ll) 



whereas the complete calculation |H| gives 

\/nl\ ~ - !)• The estimate (|B-11| ) 

reproduces the result of the complete calculation up to an order of magnitude since 

K - i|jv e = 0(1). 
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